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ABSTRACT 

We present improved calculations for the Stark broadening of hydrogen lines 
in dense plasmas typical of white dwarf atmospheres. Our new model is based 
on the unified theory of Stark broadening from Vidal, Cooper, & Smith. For 
the first time, we account for the non-ideal effects in a consistent way directly 
inside the line profile calculations. The Hummer & Mihalas theory is used to 
describe the non-ideal effects due to perturbations on the absorber from protons 
and electrons. We use a truncation of the electric microfield distribution in the 
quasi-static proton broadening to take into account the fact that high electric 
microfields dissociate the upper state of a transition. This approach represents a 
significant improvement over previous calculations that relied on the use of an ad 
hoc parameter to mimic these non-ideal effects. We obtain the first model spectra 
with line profiles that are consistent with the equation of state. We revisit the 
properties of DA stars in the range 40,000 K > Tpg > 13,000 K by analyzing the 
optical spectra with our improved models. The updated atmospheric parameters 
are shown to differ substantially from those published in previous studies, with a 
mean mass shifted by +0.034 M©. We also show that these revised atmospheric 
parameters yield absolute visual magnitudes that remain in excellent agreement 
with trigonometric parallax measurements. 

Subject headings: white dwarfs — stars: atmospheres — line: profiles 



1. INTRODUCTION 

The most successful technique used to determine the effective temperatures and surface 
gravities of hydrogen-line DA white dwarfs is to compare the observed and predicted hydro- 
gen line profiles. This so-called spectroscopic technique was first applied to a large sample of 
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DA stars by lBergeron et al.l (jl992l . hereafter BSL92) in the case of the hydrogen Balmer hnes. 
In recent years, th e method has also been apphed to the s tudy of the Lyman hne profiles 



in the ultraviolet (IBarstow et al.l l2003l : IVennes et al.ll2005l ). The success of this approach 



resides in the fact that t he theoretical line profi les are very sensitive to variations of the 



atmospheric parameters (jWegner fc Schula Il98ll ) . This is illustrated in Figure [T] where the 



theoretical profiles of five Balmer lines are displayed for various values of Tes and logg. In 
comparison with the other fitting methods used before (see BSL92 for an extensive review), 
the spectroscopic method has the lowest intrinsic uncertainties, which allows for a more pre- 
cise comparison, at least in a relative sense, of the atmospheric parameters between different 
DA stars. For instance, BSL92 used the spectroscopic method to dete rmine t he shape of the 



DA ma ss distribution (see also lLiebert et al.ll2005l : iKepler et al.ll2007l ). while iBergeron et al. 



(Il995bl ) applied the same met hod to define the locat ion of the ZZ Ceti instability strip in 
a log^f — Teff diagram (see also iGianninas et al.ll2006l ). Since about 80% of the white dwarf 
population is of the DA type, the spectroscopic technique coupled with hig h signal-to-noise 



spectroscopic observations o f the Balmer lines for large samples of DA stars (IGianninas et al. 



20051 : lEisenstein et al.ll2006l ) can reveal important details about the luminosity function and 
the evolution and the formation rate of DA white dwarfs (ILiebert et al.ll2005l ). 



Even though the relative accuracy of the spectroscopic method is mostly limited by 
the quality of the observations (signal-to-noise ratio, fiux calibration, etc.), the absolute 
values of the atmospheric parameters depend highly on the level of sophistication of the 
physics included in the calculations of model atmo spheres with hyd r ogen-r ich compositions. 
One starting point would be the efficient code of IWesemael et al.l (Il980l ) that allowed for 
the computation of radiative LTE atmospheres for hot DA white dwarfs. Similar codes 
have also been developed by D. Koester and described in iFinley et al.l (119971 ). Since then, 
the mo st important advance h as certainly been the inclusion of NLTE effects — in the 
code of iHubeny fc Lang (119951 ) among others — which allowed the study of hot DA stars 



(Tpff > 40,000 K). In terms of cooler white dwarfs, the code developed by [Bergeron et al. 
(jl99ll ) made great advances and has been a referenc e ever since. It in cludes convective 



energy transport and the non-ideal equation of state of iHummer &: Mihalad (jl988l . hereafter 
HM88) to describe the atomic populations. Otherwise, it is mostly in the atomic physics, 
with new or improved opacities and partition functions, that advances have been made. 

The dominant features observed in DA spectra are the hydrogen lines. It is thus im- 
portant to have a good understanding of the atomic transitions to exploit the power of 
the spectroscopic technique to its fullest. The main source of line broadening in most DA 
stars (Teflf > 10,000 K) is due to charged particle interactions, also called Stark broadening. 
The theory that has been the most successful to describe these line profiles is the unified 
theory of Stark broadening from Vidal, Cooper, & Smith (1970, hereafter VCS; see also 
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Smith et al. 


1969; 


Vidal et al. 


1971. 


1973) 



the analysis of large samples of DA white dwarfs (BSL92) and it is the basis of most of the 
widely accepted results on the global properties of these stars. Even though the theoretical 
framework of the unified theory has b een known since the original work of VCS, it is not 



until the calculations of iLemkd (119971 ) that a complete grid of hydrogen line profiles from 
the unified theory became available, covering the full range of possible transitions, temper- 
atures, and electronic densities encountered in DA white dwarf atmospheres. Before that, 
line profiles from the unified theory were either extrapolated to high densities (T. Schoning 
and K. Butler, private commuri i cation) , or more dras tically, less accurate broadening theo- 



ries (lUnderhill &: Waddelllll959l : Edmonds et al.l 119671 ) were used for the highest lines of the 



Balmer series (He and above). At T^s below ~ 10,000 K for the Balmer lines and below 
30,000 K for the Lyman lines, other types of line broadening, mostly due to neutral parti- 
cles, mu st be included in t he model atmospheres. The modeling of the quasi-molecular line 



opacity (lAllard et al.ll2004l ) has been one of the most significant advance in this domain. 



Despite the development of these theoretical tools, the study of DA stars has suffered 
from several complications. For instance, in their preliminary analysis of 129 DA white 
dwarfs, BSL92 discovered a lack of internal consistency between the spectroscopic solutions 
obtained when an increasing number of Balmer lines were included in the fitting procedure. 
This is illustrated at the top of Figure [2] for a typical DA white dwarf where we can see that 
the solution drifts in the Tes vs logg diagram as more lines are included in the fit. This 
correlation ind icates tha t the p hysics included in the model calculations needed significant 
improvement. iBergeronl ( 19931 ) traced back the p roblem t o the neglect of non-ideal effects 
inside the Stark broadening calculations. Indeed, ISeatonI (Il990l ) had argued that non-ideal 
effects, such as those taken into account in the HM88 equation of state, have also to be 
included directly in the line profile calculations to get a coherent physical framework. How- 
ever, it was not possible at the time of the analysis of BSL92 to rework the line broadening 
calculations. Instead, the authors chose to include an ad hoc param eter inside the model 
atmosphere code to mimic the non-ideal effects in the line profiles (IBergeronl Il993l ). By 
taking twice the value of the critical electric microfield (/3crit, see § 2.3) from the HM88 the- 
ory, it was found that the internal consistency improved substantially, as can be seen in the 
middle panel of Figure [21 It should be stressed, however, that this does not imply that HM88 
underestimated the value of the critical field. This is just a quick and dirty way to simulate 
the non-ide al effects by red ucing the line wing opacity, in particular where the line wings 
overlap (see IBergeronl Il993l for further details). This ad hoc parameter has b e en used ever 



since in all model spectra of white dwarf stars (BSL92; iHubeny fc Lanzlll995l : iFinley et al. 
I997I : iKawka fc VennesI 120061 ) . 



In this paper, we solve the problem discussed above in a more elegant way, by computing 
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improved Stark broadening profiles based of tlie unified theory of VCS, but by including 
non-ideal effects direct ly in tfie line p rofile calculations following an approach similar to that 
originally proposed by ISeatonI (Il990l ). The non-ideal effects due to the proton and electron 
perturbations are taken into account using the HM88 theory. We begin in § 2 by discussing 
the two successful theories at the heart of our calculations, namely the VCS unified theory 
of Stark broadening and the HM88 non-ideal equation of state. In § 3, we describe how these 
two formalisms are combined for the first time to compute Stark broadening profiles that 
take into account non-ideal effects. We then evaluate in § 4 the implications of our improved 
line profiles on the spectroscopic analysis of DA stars. Our conclusions follow in § 5. 



2. LINE BROADENING THEORY 



The basis of model atmospheres for DA white dwarfs resides in the equation of state and 
opacity calculations of the hydrogen gas. One of the most successful equations of state used 
to compute the populations of the different accessible states (bound states of hydrogen, H~, 
molecular H2, and H.^, as well as proton and electron populations) is the HM88 non-ideal 
equation of state described in more detail in section § 2.2. With these populations in hand 
for all states of the gas, we can proceed to calculate the gas opacity, which is generally split 
into true absorption and scattering processes. In particular, the bound-bound opacity (or 
line opacity) between one initial state i and one final state j is written as 



ire 



ITLpC 



where Ni is the population of the initial state and fij is the oscillator strength of the tran- 
sition. These two factors define the amplitude of the transition while 0(i^) is the spectral 
broadening profile normalized to unity. In this paper, we are mostly interested in the broad- 
ening profile for the hydrogen lines. This opacity is the key ingredient of the spectroscopic 
method that relies on a detailed comparison of observed and predicted line profiles. 

Generally speaking, the line profile is a convolution of different broadening mechanisms. 
The first source of broadening is due to the interaction between the absorber and the charged 
particles in the plasma. This process is referred to as Stark broadening and is discussed in 
more detail in section § 2.1. Stark broadening represents the dominant source of broadening 
for most transitions in DA w hite dwarfs. Accounting only for this source of br oadening, the 
line profile can be expressed (lUnderhill fc Waddelllll959l : Edmonds et al.l 119671 ) as 
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\2 

(l)(u)du = —S*(a)du , (2) 
cFo 

where a = AX/Fq with A A measured from the center of the hne. Fq is defined as the electric 
field at the mean distance between the plasma ions 



Fo = Zpe(47rArp/3)2/3 = 1.25 x IQ-^ZpN^^^ , (3) 

where Zp and Np are the charge and density of perturbers, respectively. The profile S*{a) 
is a convolution of a pure Stark profile S{a) and a Voigt function H{a,v) that takes into 
account thermal and natural broadening 



S*(a) 



S a + 



\ov\ H{a,v] 



dv 



(4) 



where a = r/47rAz/£), with F representing the natural broadening half- width and Ai>£, the 
Doppler width. The second type of broadening important in DA white dwarfs is the in- 
teraction of the absorber with neutral hydrogen. It is especially important at low effective 
temperatures for transitions involving an upper level with a low pr incipal quantum number . 
In the case of Balmer lines, w e must include resonaii ce broadening (lAli fc Griemlll965l . Il966l ) 
and nonresonant broadening (IHammond et al.lll99lh in cool white dwarfs with Tcs < 10,000 
K. Since this source of broadening has a Lorentzian profile, the broadening parameter can 
be added directly to that of the natural broadening. For the lower Lyman lines (La, LjS, and 
L7), close range collisions of the absorber with hydrogen atoms and protons cause the ap- 
pe arance of important sate l lites in the wings o f the lines that are visible up to Tgg ~ 30, 000 
K f AUard fc Kielkopl 1982 ; Allard et al.lboO-j ). This opacity also affects the thermodynamic 
structure of the atmosphere. Furthermore, at v ery low effective temperatur es, H-H2 collisions 
become the main source of broadening for La ( iKowalski &: Saumonll2006l ). 



2.1. Stark Broadening 

The Stark effect is defined as the shifting — or splitting — of spectral lines under the 
action of an electric field. In the following section, we assume that we have in the atmo- 
spheric plasma a local electric microfield due to protons that is constant with time. This 
static approximation is well justified since the characteristic time for the fluctuation of the 
local microfield due to a change in the proton distribution is much large r than the character- 



istic time of absorption and emission processes that we are studying (IStehle &: Jacquemot 
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19931 ). We also assume that the microfield is constant in space, a good approximation if the 



perturbers are far from the area occupied by the bound electron orbits. The regime where 
this approximation fails is discussed further in § 2.3. The electric microfield, assumed to be 
along the z-axis, interacts with the dipolar moment of the atom and leads to a perturbation 
in the Hamiltonian of the form 



H, = eR-F = eFz . 



(5) 



It is common to write the amplitude of the microfield in terms of the unitless parameter 
j3 such that 



(6) 



where Fq is the characteristic field defined by equation (j3]). The amplitude of the microfield 
can be expressed using a probability distribution P{P). One such distribution that can be 
computed analytically is the Holtsmark distribution, which takes into account the inter- 
actions (vecto rial sura ) betw een the different perturbers. We use here the more physical 
distribution of iHooperl (119681 ) that takes into account the Debye screening effect, and which 
corresponds also to the formulation included in the work of VCS. H M88 originally u s ed the 
Holtsmark distribution, but for consistency we use the approach of iNayfonov et al.l (119991 ) 
for the HM88 equation of state by replacing the Holtsmark distribution with the Hooper 
distribution. We note that this improved version of the HM88 equation of state is already 
included in our model atmosphere code, although it has never been properly documented. 

We note that the perturbation Hg does not commute with the angular momentum 
since it has a privileged direction in space. The effect of breaking this symmetry is to lift 
the degeneracy in energy of the hydrogen atomic levels. There are no exact solutions to 
this problem and we generally use quantum mechanical perturbation theory to determine 
corrections to the energy and wave functions. The first order corrections to the energy are 
linear in relation to the electric microfield that is applied, and represent the basis of the linear 
Stark effect. The linear corrections are found by diagonalizing the perturbation operator in 
the degenerate subspaces. For a hydrogen atom, we can write this expression analytically. 
First of all, knowing that the operator z conserves the azimuthal q uantum number m, we 
can diagonalize the matrix (\Effc| if^ (ICondon fc Shortleyl 119351 ) by making use of the 
quantum numbers {n,m,q), where q = qi — q2 with gi + ^2 = n — \m\ — 1, and n is the 
principal quantum number. The shifts in energy are then given by 



(7) 
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When the electric microfield becomes very high, or equivalently, when the Stark sphtting 
is important enough that the Stark components of two levels with a different principal 
quantum number are crossing, the linear approximation is no longer valid. The second order 
corrections to the energy become quadratic in relation to the microfield, and we find 

A 77.2 ^2z?2 I (^n',m',g'| Z |\&n,-m,g) P /o\ 
^^n,m,9 = e ^ 2^ • (8) 

n'^n n,m,q ^n',m',q' 

The problem with this perturbation approach to the Stark effect is that the wave func 



tions, using any order of the perturbation theory, are not necessarily normalizable (IFriedrich 



20061 ). This is because the perturbing potential goes to —00 when z goes to —00. In other 
words, states that were bound without a microfield are now only metastable states, and 
there is a finite probability that the atom will be ionized when the microfield has a suffi- 
ciently high amplitude. Classically, the sum of the electric potentials of the absorber and 
the nearby protons only allows bound states in the local potential minima up to a certain 
energy, called the saddle point. Part of this problem is because we have used a spatially 
uniform microfield (eq. E]). In reality, the microfield cannot be uniform beyond a distance 
of about the interparticle separation, and thus under normal conditions, most of the states 
will be bound. However, for sufficiently high microfields, some excited states will become 
unbound and this is where resides the problem discussed further in the next sections. 



2.2. The Unified Theory of Stark Broadening 



In addition to static protons, the absorber also interacts with free electrons. These 
particles being much faster than the protons, a coUisional approach is generally used to 
describe this interaction. Since the collisions are rapid, the interaction is non-adiabatic and 
a proper quantum mechanical treatment of the collisions accounting for the internal structure 
of the hydrogen atom is required. Various theories of Stark broadening for hydrogen lines 



take into account both the electron and proton interactions (IVidal et al.lll97ll : ISeatonlll990 



Stehle &: Jacquemotlll993l ). Here we summarize the principal aspects of the unified theory 
of VCS, which even today stands as the most accurate theory for the conditions encountered 
in white dwarf atmospheres. 

The unified theory of VCS uses the quasi-static proton broadening approximation with 
the distribution of electric microfields described in § 2.1, assumed to be constant in space 
and time. The complete Stark profile is then defined as 
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S{a) = / P{P)I{a,P)dp , (9) 
Jo 

where I{a,f3) is the electronic broadening profile. Therefore, the problem is reduced to the 
calculation of electronic broadening profiles for a fixed microfield of amplitude /3. The full 
broadening profile is then the average of the electronic profiles over all possible microfields 
weighted by the probability distribution -P(/3). The unified theory is a non-adiabatic quan- 
tum mechanical theory that uses the classical path approximation for the electron-absorber 
interactions, implying that the wave functions are well separable. The interaction potential 
considers only the first dipole term. The name of the unified theory comes from the fact 
that the electronic profile is valid at all detunings from the line center. In the asymptotic 
limit of the line core, the impact approximation is recovered. Furthermore, in the far wings 
of the lines, the one-electron theory is recovered and the electrons and protons have similar 
contributions to the broadening. 

Now we want to compute the electronic profile for an initial Stark state to a final state 
n'^, where this abbreviated notation is used to designate any of the [n, m, q) states. Mostly to 
reduce the computing time and to simplify the inversion of the matrix involved (see below), 
the unified theory uses the no-quenching approximation, which means that there are only 
collision-induced transitions within the same principal level. The states after the collisions 
will be denoted by and nj, for the initial and final state, respectively. The calculations 
can then be performed by considering only two-level transitions. This last approximation 
is valid only if the difference in energy between the highest Stark state from level n' and 
the lowest state of level n' + 1 is large, or in other words, when the different lines are well 
separated (i.e., when the line wings do not overlap). Only then is the regime of the linear 
Stark effect consistent with the basic assumptions of the unified theory. 

The electronic broadening profile I{a,(3) is given by the sum over all states 



IT 



pina)lm {ual d\n'^) {n'f^l d\nb) {nb\ {n'f^l )C ^{a, (3) \n'^) \na) 



(10) 



where d is the dipole operator, and p is the probability that the atom is in the initial state 
(Boltzmann factor). The operator / C takes into accou nt the linear Stark effect as well as the 



free electron-atom interactions (see lVidal et al.lll970l . eq. XII. 2 for more details). Since the 



operator K. is generally not diagonal in the basis of the (n, m, q) states, it induces a coupling 
between the various Stark components. 
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2.3. Non-Ideal Equation of State 

One important advance tliat lias provided a better interpretation of the h ydrogen lines in 



DA w hite dwarfs is a realistic modeling of non-ideal effects at high densities. iBergeron et al. 



( 119911 ) were the first to include the occupation probability formalism of HM88 to deter- 
mine the populations of the bound states of hydrogen in white dwarf atmospheres. This 
probabilistic approach considers the perturbations on each atom by charged and neutral 
particles. One advantage of this statistical interpretation is that there are no discontinuities 
in the populations and the opacities when the temperature and the pressure vary, in contrast 
with other formalisms that simply predict the last bound atomic level under given physical 
conditions. 

Briefly, each atomic level n has a probability Wn of being bound and a probability 1 — w„ 
of being dissociated due to perturbations from other particles in the plasma. The partition 
function of a given species is then written as 

Z = ^Wngnew(y-^) , (11) 

n 

where Xn and Qn are respectively the excitation energy and multiplicity of the level. The 
occupation probability w„ is defined as 

= exp ( 1 , (12) 

where / is the free energy of the non-ideal interaction. The different interactions being 
statistically independent, the total occupation probability can then be calculated simply as 
the product of the contributions from neutral and charged particles. 

The interaction with neutral particles is treated within a hard sphere model. The 
configurational free energy is derived from the second virial coefficient in the van de Waals 
equation of state (excluded volume correctio n). The exact origin of this excluded volume 



term is the hard sphere equation of state of ICarnahan fc Starlingl (119691 ). The interaction 



with charged particles is the one closely related to the Stark effect. This interaction suggests 
that the electric microfields, fluctuating on long time scales due to changes in the spatial 
distribution of protons, can destabilize a bound state and cause its dissolution into the 
continuum. Indeed, we already discussed in § 2.1 that the application of an intense electric 
field on a bound atom allow for the ionization of the atom. HM88 therefore suggest an 
occupation probability for proton perturbations of the form 
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w„ (charged) 



P(/3) dp, 



(13) 



where P{f3) is the probabihty distribution of the electric microfields introduced in § 2.1. This 
occupation probabihty imphes that all microfields larger than the critical field /3crit will ionize 
the electrons in level n (it is implicit here that the critical field depends on the atomic level 
considered). The difficulty with this formulation is to find the value of the criti cal microfield 



for which slightly differen t expressions can be found in the lit erature (HM8 8; ISeatorull990 



Stehle fc Jacquemotlll993l ). The simplest formulation is that of ISeatonI (jl990l ) which consists 
in taking the critical microfield as the point where the energy of the highest Stark state 
for a given level n crosses the energy of the lowest Stark state of the next level n + 1. For 
hydrogen, and using the linear Stark effect, we find 



_ 2n + l e 
~ 6n^{n + 1)2 alF, 

The assumption that the crossing of two atomic levels with different principal quantum 
numbers leads to the dissolution of the lower level is difficult to prove and it is based in part 
on laboratory experiments and theoretical considerations (see HM88 for a more detailed 
discussion). Let us consider one bound electron in one of its many Stark states of the n- 
th level. When the electric microfield, which fiuctuates in time, gets to a value such that 
there is a crossmj^ between this level and the lowest Stark state of the n + 1 level, there 
is a significant probability that when the microfield goes down, there will be a transition 
n —>■ n + 1. The incessant fiuctuations of the microfield imply that these transitions will 
continue until the electron becomes unbound. However, the bound electron will in fact be 
in a rather homogeneous superposition of all its accessible Stark states due to electronic 
collision^ on time scales much shorter than the fiuctuation time for the microfields, and also 
due to fiuctuations in the direction of the microfield. That is, as soon as the first Stark state 
from a level n crosses another one from the level n + 1 , the electron is allowed to cascade to 
the continuum. 

The reader might have noticed that equation fll4p is based on the linear Stark effect 
although this is not necessarily a good approximation when levels are crossing. This is 
why HM88 considered a formulation of the critical field that takes into account corrections 



^Strictly speaking, direct degenerate crossings are avoided due to fine structures and non-uniformity in 
the microfields. 

^See equation [ini the non-diagonal operator JC introduces this coupling. 



(14) 
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of higher order from the perturbation theory as well as results from experiment^. They 
find that for levels with n > 3, the linear theory remains valid. For lower levels, however, 
the critical field is large enough that it must have been created by a single proton very 
close to the absorber. In this case, the approximation of a uniform electric field in space 
fails and we must consider the full potential curve of the system to find the critical field 



( Stehle &: Jacaueniotlll993h . and also to compute the broadening profiles of the far line wings 
(lAUard &: Kielkopj|l982l ). Exact calculations show that no crossing is possible between levels 
n = 1, 2, and 3, which implies that La is not affected by non-ideal effects due to proton 
perturbations. Also, the crossing between n = 3 and 4 is close to the classical saddle point 
value, which means that the linear Stark theory is in error by ~ 14%. Therefore, HM88 
added a smoothing factor to equation (|T^ so that the exact result is recovered for n = 3. 



For levels with a high value of n, HM88 obtain the same results as ISeatonI (Il990l ). For 
the Balmer lines, these different values of /5crit are not an issue since the non-ideal effects 
discussed above are never important for Ha in DA white dwarfs. This is not the case for 
La and L/3, however, for which non-ideal effects should be treated with caution, or simply 
neglected (for instance in TLUSTY and in our code discussed in § 4). 

We have seen in § 2.2 that the interaction of hydrogen atoms with rapid electrons 
was better described with a collisional approach since the quasi-static theory is generally 
not valid in this case. These perturbations are also considered in the HM88 model. We 
must keep in mind that because the plasma is in thermal equilibrium, the equation of state 
already accounts for the inelastic electronic collisions that cause transitions and ionization. 
However, we expect that if the time between collisions is small, the energy of a bound level 
will be left undefined by a certain amount given by the uncertainty principle. If this energy 
uncertainty is of the order of the ionization potential, we expect the state to become unbound 
after the collision. As stated in HM88, this contribution to the occupation probability is 
always much lower than the proton contributions and can be neglected in DA atmospheres. 
However, the electronic perturbations should not be neglected when calculating the line 
profiles, as discussed in § 3. We must also remember that the electrons contribute indirectly 
to the occupation probability due to proton perturbations by mixing the Stark states. The 
occupation probability for the bound states of hydrogen are displayed in Figure [3] for typical 
conditions encountered at the photosphere of white dwarf stars. Non-ideal effects are shown 
to be extremely important except for the very lowest transitions. The solution of the HM88 
equation of state yields the populations for all states. 

Finally, the opacity calculation must also be modified to take into account the occupation 



•^Remember that even if we have the exact critical field, ionization docs not necessarily occur exactly at 
that point and this is why experimental data are also required. 
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probability formalism, as described by lDappen et al.l (119871 ). The bound-bound (line) opacity 
for a transition between levels i and j (equation [1]) is multiplied by Wj/wi to take into account 
the fact that these levels can be dissolved when a photon is absorbed. We thus obtain 



Tie w 



nif. w. 



^ fij<p{v) dv 



(15) 



When the absorption is from a bound state to a final state that has been dissolved by non- 
ideal interactions (with a probability 1 — Wj/wi), we obtain instead a bound-free opacity 
since this process is equivalent to the ionization of the atom. We can directly extrapolate, 
in this case, the bou nd-free cross section at frequencies below the usual cutoff frequency Uc- 
Dappen et al.l (119871 ) have treated this opacity by considering the absorption from a level i 
of a photon of energy hu that yields a transition to a fictitious upper level n* given by 



n 



The hydrogen bound-free opacity for v < 
can then be written as 



hv 
7 



-1/2 



(16) 



also called the pseudo-continuum opacity 



Ki{u)du = A^j I 1 — 



Wi 



■ du 



;i7) 



In Figure m we show the individual contributions of the line and pseudo-continuum opacities 
at the photosphere of a 20,000 K DA white dwarf in the spectral region of the Balmer lines. 
We can see how the pseudo-continuum opacity extends to wavelengths longward of the 
Balmer jump, potentially affecting the opacity between the lines, and in particular the high 
Balmer lines. 



3. STARK BROADENING PROFILES INCLUDING NON-IDEAL 

EFFECTS 



The VCS line profiles (§ 2.2) and the HM88 equation of state 2.3) have long been in- 
cluded in modern white dwarf model atmospheres (BSL92; Hubeny fc Lanj[l995 : Finley et al. 
19971 : iKawka fc VennesI l2006l ) . The HM88 theory has been used, in particular, to calculate 
populations as well as line and pseudo-continuum opacities. However, all model spectra 
currently availab l e and calculated within the HM88 framework suffer from a major incon- 
sistency. ISeatonI (jl990l ) was the first to point out that not only the line strengths will be 



- 13 - 



affected by non-ideal effects, but tlie shape of tlie line profiles as well. The Stark broadening 
profiles S{a) from VCS, for instance, do not take into account such non-ideal effects. It is 
possible to describe qualitatively the behavior of non-ideal line profiles. The high electric 
microfields, normally responsible for the absorption in the far wings of the line profiles, will 
not contribute as much to the line opacity in the non-ideal case since these fields are mostly 
responsible for dissolving the atomic levels. Consequently, the opacity in the line wings will 
be significantly reduced with respect to the ideal case. But since by definition the profiles are 
normalized to unity, the non-ideal line profiles will have deeper cores and appear narrower 
in comparison with the ideal case. 



Instead of including the non-ideal effects discussed above, iBergeronl (Il993l ) proposed 
a solution that could be easily implemented in model atmosphere calculations. As men- 
tioned in the Introduction, this is the solution adopted by BSL92 in their analysis of a large 
spectroscopic sample of DA stars. The authors found in their preliminary analysis that the 
atmospheric parameters varied significantly when an increasing number of Balmer lines was 
included in the fitting procedu re. Furthermore, the m ean mass of their sample w as uncom- 
fortab ly low ((M) ~ 0.53 Mq] iBergeron et al.lll990al ). The solution proposed by [Bergeron 
(I1993I ) was to parameterize the value of the critical field in equation (fHll . The effect of 
varying the value of Pcrit is, among other things, to change the relative importance of the 
line wing and pseudo-continuum opacities. By increasing arbitrarily the value of the criti- 
cal field, one mimics the non-ideal effects by reducing the line wing opacity. One problem 
with this approach, of course, is that the determination of the multiplicative factor in the 
expression for /3crit has absolutely no physical basis. This factor can only be determined 
from an empirical analysis of t he internal cons istency of the fitting procedure, similar to 
that shown in Figure [2] (see also IBergeronl Il993l ). Furthermore, there is no reason to expect 
this ad hoc parameter to be the same for all hydrogen levels (i.e., Balmer or Lyman series) 
or for different atmospheric parameters. Another side effect of this approach is to artificially 
change the atomic populations in the HM88 equation of state, which in turn could have an 
unexpected impact on other parts of the model calculations. It should be clear by now that 
a more physical approach is seriously required, especially given the high quality spectra that 
have become available in recent years, either in the visible or the ultraviolet. This ad hoc 
parameter represents an important hurdle towards our understanding of the global properties 
of white dwarfs. 



A coherent way to comb ine the Stark broadening profiles and the HM88 equation of 
state is described at length in ISeatonI (Il990l ). Seaton proposed to replace equation by 
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S{a) 



Z^crit 



P{(3)I{a,(3)d(3 



/3crit 



P{l3)dj3 



where only electric microfields j3 with an amplitude inferior to the critical field jScnt for the 
upper level of the transition now contribute to the broadening of the lines by the protons. 
The higher fields contribute instead to the pseudo-continuum opacity. The denominator, 
which turns out to be the occupation probability due to proton perturbations, allows for the 
renormalization of the line profile^. 



SeatonI (119901 ) was the first to calculate Stark broadening profiles taking into account 



non-ideal effects. However, his calculations made for the Opacity Project rely on an ap- 
proximate ele ctron broadening the o ry th at is inappropriate in the context of white dwarf 
atmospheres. IStehle &: JacquemotI (1l993l ) also discuss the implementation of non-ideal ef- 
fects inside the line profile calculations using an alternative electron broadening theory. This 
theory also lifts in part the static approximation for the protons and therefore gives better 
results in the center of the profiles, but this has no effect on the spectrum at the high densities 
found in white dwarf atmospheres. In the line wings, the results are equivalent to the VCS 
formulation. We must men t ion, however, that all pu bhshed tables and subsequent papers 
for this theory (jStehlelll994J : IStehle &: Hutcheonlll999l ) do not include non-ideal effects. 



In this work, we use the Seaton approach (equation [18]) to include non-ideal effects 
due to proton perturbations inside the unified theory of Stark broadening from VCS. We 
also include non-ideal correct ions due to e lectronic perturbations according to the HM88 
theory. This was neglected by lSeatonI ([19901), likely because it was estimated that the effects 
would be negligible in the context studied. As discussed in § 2.3, the collisions between free 
electrons and the absorber leave a bound level in a state of indefinite energy. Therefore, the 
HM88 model states that if the uncertainty on the energy becomes equal to the ionization 
energy of the level, the electron has a good probability of becoming unbound. Here we have 
at hand the electronic broadening profiles given by equation ([TU]) . and their wings generally 
extend to infinity in terms of the detuning Az/ (measured from the line center). However, 
according to the HM88 theory, it is obvious that the detuning induced by the electronic 
collisions cannot extend much beyond (i-e. the ionization threshold frequency for this 
level) or there would be line opacity shortward of the Balmer or Lyman jumps. An abrupt 



"^In the definition of Seaton ( 1990l ). the hne profiles are not renormahzed but instead the factor Wj is 
omitted in equation [T5l 
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cutoff is not desirable since this could cause discontinuities in the spectra. We adopt here an 
exponential cutoff so that the electronic broadening profiles are reduced by a factor of 
at Au = z/c — z/q, where uq is the central frequency of the transition, and then we renormalize 
the total broadening profiles. Since the exponential factor is independent of f3, we can easily 
represent this correction in the form 

</>'(^) = • (19) 



This has a small but non-negligible effect on the line wings, especially for the higher members 
of the series. Both Tes and logg are increased by ~ 0.5 % when this correction is introduced 
in our calculations. Using alternative types of cutoff, we estimate that the uncertainties on 
this part of the theory can be as high as half of the shifts found here. As mentioned in 
HM88, a more physical approach is certainly much needed but such a theory still does not 
exist. 

Figure [5] presents the Stark broadening profiles S{a) for the Balmer lines H7 and H8 at 
typical photospheric conditions of DA stars. Al ong with the results from this work, we show 



the or iginal calculations of VCS (jLemkd 119971 ) and the approximate profiles from ISeaton 



(jl990l ). We note that our profiles including non-ideal effects due to proton and electron 
perturbations are very different from the original ideal gas calculations of VCS, even though 
we are using the same broadening theory. Our profiles have significantly less opacity in the 
line wings, and higher lines of the series have much sharper profiles, as anticipated. To have 
a better idea of the relative importance of the two improvements discussed above, we have 
estimated that the electronic corrections account for 12% and 16% of the reduction of the 
opacity at 100 A from the line center of H7 and H8, respectively Finally, we show that our 



profiles are not in agreement with the calculations of ISeatod (Il990l ). which include non-ideal 
effects due to proton perturbations but use only an approximate treatment for the electronic 
broadening. These results confirm that this last approximation is not suitable for typical 
conditions encountered in DA white dwarf atmospheres. The comparison of our Lyman line 
profiles with the calculations of VCS or Seaton is similar to those displayed in Figure [51 

We have demonstrated in the previous sections that our approach to Stark broadening 
provides significant physical improvements over previously available calculations. When 
included in a model atmosphere, the equation of state and the opacity cal culations becom e 



entirely consistent without the need of the ad hoc parameter introduced by iBergerod (119931 ). 
By including non-ideal effects, we have shown that the line profiles from the unified theory 
are dramatically modified in typical white dwarf atmospheres. These improved profiles now 
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need to be validated using spectroscopic data. 



3.1. Comparison with Laboratory Experiments 



The preferred way to validate Stark broadening theories f or hydrogen has b een the 
comparison with laboratory plasma emissivity measurements from lWiese et al.l (119721 ). These 
pure hydrogen plasma arc experiments were performed at high densities and temperatures, 
comparable to the photospheric conditions of cool white dwarfs. The main advantage of this 
experiment is that it is not under the constraint of radiative equilibrium, unlike a stellar 
atmosphere. In other words, the emergent flux at one wavelength is not affected by what 
occurs at other wavelengths. The emissivity can then be simply written as 



(20) 



where k{X) represents the total monochromatic opacity. Another advantage, at least from 
a theoretical point of view, is that it is easier to control the experiments and the results 
are independent of various sources of uncertainty intrinsic to DA model atmospheres (i.e., 
convection, contamination from heavy elements, etc.). We use the data for the lowest and 
highest density experiments (the most extreme positions on the p lasma arc) frorn Wie s e et al., 
which were conapared to broadening theories in v arious studies (jPappen et al.l 119871 : ISeaton 



1990l : iBergeronl Il993l : IStehle fc JacquemotI Il993l ). Wiese et al. estimate the LTE plasma 
parameters for these two experiments at logT = 4.00, logiVg = 16.26 and logT = 4.12, 
logiVe = 16.97, respectively. However, there is no simple way to measure the plasma state 
parameters and it is not clear whether LTE is reached or not. The authors used the total 
intensity from two lines, together with the intensity at two continuum points, and fitted 
the data with an approximate plasma model. Since then, slightly different parameters have 
been used to compare with broadening theories. Instead of choosing approximate values 
for the parameters like in previous analyses, we performed a fit to the full data sets, 
using the same input physics as for our white dwarf models. The results with the original 
VCS calculations and our improved line profiles are displayed in Figures and [7] in terms of 
absolute predicted fluxes (i.e., without any renormalization) . 

For the lower density experiment of Figure [6l our improved profiles provide a much 
better fit to the laboratory data, in particular in the regions of the high Balmer lines where 
the wings overlap. The higher density experiment of Figure [7] is more problematic. We find 
that there is a partial T — N^. degeneracy in the diagram; it is indeed possible to find 
many acceptable solutions from visual inspection by increasing both T and N^- However, 
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to provide the best overall fit with our new profiles, we have to increase the plasma state 
parameters significantly compared to the Wiese et al. values. We see that our profiles provide 
the best fit for the higher lines, although the red wings of H/5 and H7 are predicted a bit 
too weak. All in all, we conclude that the experiments of Wiese et al. do not provide such 
a stringent constraint on the broadening theories because of the large range of acceptable 
plasma state parameters, and also because of potential departures from LTE. 



3.2. Further Theoretical Improvements 



While our work represents a significant advance over previous calculations, second order 
effects in the theory of Stark broadening could still come into play. Such second order 
effects have been stu died extensively i n the literature in the case of high den sity hydrogen 
plasma experiments (ILee fc Okslll998l : IStehle et al.l 120001 : iDemura et al.ll2008l ). The sources 
of these second order effects are plentiful but they all have the same impact on the line 
profiles: they produce asymmetries and they shift the central wavelengths, the latter being 
unimportant for the spectroscopic te chnique used to meas ure the atmospheric parameters 
of white dwarfs. The recent study of iDemura et al.l (120081 ) reveals that the various second 
order effects compete and interact together in a complex way and that they should all be 
included in the models simultaneously. The most important effects are due to quadrupole 
proton-absorber interactions and to the quadratic Stark effect, but second order effects due 
to quadrupole electron-absorber interactions are noticeable as well. The lifting of the no- 
quenchin g approximation c ould also be considered a second order effect, although it was 
shown by lLee fc Okd (119981 ) to have an impact mostly on the line cores. Furthermore, some 
of these quenching effects have already been accounted for implicitly by including here the 
HM88 non-ideal effects directly inside the line profile calculations. Finally, even if we include 
better physics for the atomic transitions, we still have to rely on the HM88 equation of state, 
which is, or may become, the main source of uncertainty. 

Our analysis above of the Wiese et al. experiments, with our improved profiles has al- 
ready revealed the existence of line profile asymmetries. For the lower density experiment 
(see Fig. [6]), it is only a mild effect in the far wings, however. It is not clear if such second 
order effects would be significant in a typical DA white dwarf at logg ~ 8, for which the 
electronic density at the photosphere is roughly halfway between the two Wiese et al. ex- 
periments. Therefore, it is probably premature at this stage to include second order effects 
in our calculations, and it is not clear whether such changes would affect our line profiles 
significantly. We come back to this point in the analysis of the PG spectroscopic sample 
discussed below. 



- 18 - 



APPLICATION TO WHITE DWARF ATMOSPHERES 



4.1. Model Spectra 

We now discuss the astrophysical implications of our improved line profiles described 
in § 3 on the modeling of DA w hite dwarfs. Our model atna osphere code is based on the 
program originally developed by iBergeron et al.l (jl99ll . Il995al ) and references therein. The 
main difference is that we have updated several sources of opacity and partition functions 
(all within the occupation probability formalism). For completeness, we provide in Table 1 
the complete list of opacity sour ces included in our code. W e have not included the La line 
broadening due to H2 collisions (IKowalski &: Saumonll2006l ) since these calculations are still 
not available, although this has no impact over the range of effective temperature considered 
here. 

We restrain our analysis to effective temperatures above ~ 13, 000 K to avoid additional 
uncertainties related to convective energy transport. Indeed, the atmospheric structure of 
DA white dwarfs below this temperature depends sensitively on the assumed convective 
efficiency. The most commonly used approach to include convection in model atmosphere 
calculations, the mixing-length theory, is at best a very crude approximation, and the convec- 
tive efficiency niust be parameterized by carefully adjusting the value of the mixing length 
(IBergeron et al.l Il995bl ). It is also believed that large amounts of helium can be brought 



to the surface by convection whi le remaining spectroscopically invisible (IBergeron et al. 



1990bl : iTremblay fc Bergeron! |2008| ). Hence, it is not even clear whether cool DA stars have 
hydrogen-rich atmospheres. Another reason for restricting our analysis to hotter stars is 
that neutral line broadening becomes important for the Balmer lines below Tes ~ 10, 000 K, 
preventing us from performing a direct comparison of the observed and predicted Stark 
broadened absorption lines. We also set an upper limit of T^s = 40, 000 K because NLTE 
effects become important above this temperature. Also, absorption from heavier elements, 
particularly in the UV, is likely to affect the atmospheric structure of hot white dwarfs. 

We thus computed two grids of model atmospheres, one with the VCS profiles and 
the other with our improved profiles (both calculated with /?crit x 1), covering the range of 
Teff = 12, 000 K to 45,000 K (with steps of 500 K up to 15,000 K, 1000 K up to 18,000 K, 
2000 K up to 30,000 K, and 5000 K above) and of log g = 6.5 to 9.5 (by steps of 0.5 dex with 
additional models at 7.75 and 8.25). To compare the model spectra calculated from both 
sets of line profiles, we simply fit the VCS spectra with our improved spectra. The results 
of this exercise are illustrated in Figure [HI In general, our new models yield systematically 
higher effective temperatures (by 1000 — 2000 K) and surface gravities (by ~ 0.1 dex) over 
the entire range of atmospheric parameters considered here. One exception is a region at 
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low temperatures and high surface gravities where the trend in temperature is reversed. The 
strip, inchned in the diagram near T^s ~ 15, 000 K, where this separation occ urs corresponds 



to the region where the strength of the Balmer hues reaches its maximum (IBergeron et al. 



1995bl ). Our new models predict Balmer lines that are stronger, hence the temperatures in 



this region are pushed towards lower values. 



4.2. Analysis of the DA Stars in the PG Sample 

In this section, we measure the implications of our improved Stark profiles on the anal- 
ysis of Balmer line observations of DA stars by considering the Palomar-Green sample of 



Liebert et al.l (120051 ). This sample of 348 DA stars has been studied in great detail, and the 
range of effective temperatures for these objects corresponds very well to that considered 
in our analysis. About 250 white dwarfs fall in the appropriate range of Teg, depending on 
which model grid is used. The data set and the fitting pro cedure of the optical spectra are 



identical to those described at length in lLiebert et al.l (120051 . and references therein). Briefly, 



we flrst normalize the observed and model flux from each line to a continuum set to unity at 
a flxed distance from the line center. The observed proflles are then compared with the pre- 
dicted proflles, convolved with a Gaussian instrumental proflle. The atmospheric parameters 
are then obtained using the non-linear least-squared method of Levenberg-Marquardt, fltting 
simultaneously flve lines {B.(3 to H8). In some cases where the spectrum is contaminated by 
a M dwarf companion, one or two lines are excluded from the flt. 

As discussed in the Introduction, model spectra calculated with standard Stark broad- 
ening proflles yield inconsistent atmospheric parameters when different lines are included in 
the fltting procedure. As shown in the bottom panel of Figure [21 our improved Stark proflles 
provide an even better internal consistency than the previous calculations displayed in the 
two upper panels. To better quantify this internal consistency, we perform for each star in 
the PG sample the same exercise as that shown in Figure [2] using the VCS proflles and our 
improved proflles. We then compute for the 250 stars in our sample the average absolute 
deviations in T^g and logg between the solutions obtained with a different number of lines 
included in the flt. The results of this exercise are presented as fllled circles in Figure [9] to- 
gether with the mean uncertainties from the fltting procedure as a reference point (the open 
circles will be discussed in § 4.3). We can see that the deviations with our improved proflles 
have been signiflcantly reduced by a factor of ~ 1.6 in T^g and ~ 1.8 in logg. Furthermore, 
these deviations now both lie within the mean uncertainties of the fltting procedure, a result 
which is extremely reassuring. 



The quality of the flts is also an important aspect of the comparison between model 
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grids. We compare in Figure [TOl our best fit to WD 0205+250 (same star as in Fig. [2]) using 
the VCS profiles and our improved Stark profiles. This comparison reveals that while the 
atmospheric parameters are significantly different, the quality of the fits is similar. We thus 
conclude that the quality of the fits cannot help to discriminate between both sets of Stark 
profiles. 

We can also explore the global properties of the PG sample with our two model grids. We 
first convert the log g values into mass using evolutio nary models appropriate for white dwarfs 
with thick hydrogen layers (see lLiebert et al.ll2005l for details). Our results are presented in 
Figure [11] in a mass versus effective temperature diagram. As expected from the previous 
comparison displayed in Figure [TOl our improved line profiles yield significantly larger masses 
and higher effective temperatures. We note that the results displayed in the bottom panel 
are qualitativel y similar to the early s pectroscopic determination of the mass distribution 
of DA stars by [Bergeron et al.l (jl990al ). which was based on the VCS profiles without any 
modification of the critical field. In this case, the mass distribution has a mean value near 
0.53 M0, which is uncomfortably low compared to what is expected from earlier phases of 
stellar evolution (see the discussion in [Bergeron et al.lll990al ). 



Finally, we have also investigated whether second order effects in the Stark broadening 
theory, such as those discussed in § 3.2, could be detected in our analysis. In particular, if 
such effects are present, we would expect to observe differences in the fits of the blue and red 
wings of the Balmer lines. We have thus reanalyzed the PG sample by including in the fit 
(1) only the red wing of the Balmer lines and (2) only the lines cores (this is accomplished 
by fitting only half of the wavelength range we normally use for each line). We find that 
the atmospheric parameters obtained in this manner are entirely consistent, and we thus 
conclude that second order effects can be safely neglected in DA white dwarfs. 



4.3. A Reappraisal of Previous Studies of DA White Dvi^arfs 

The spectroscopic determination of the mass distribution by BSL92 was based on model 



spectra that include the solution proposed by [Bergeron! (1 1993! ) to mimic the non-ideal effects 



namely by taking twice the value of the critical field (/3crit x 2) in the HM88 formalism. This 
ad hoc procedure had the effect of artificially reducing the pseudo-continuum opacity, and 
thus the opacity in the wings of the high Balmer lines. The internal consistency of the 
solutions obtained from different lines was consequently improved, as can be judged from 
the results displayed in Figures [2] and [9[ (open circles). Since this ad hoc solution has been 
adopted by BSL92 and in all white dwarf models used in the literature, it is important to 
evaluate the differences between our revised atmospheric parameters and those published 
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in the literature based on these approximate model spectra. For this purpose , we have also 
calculated another model atmosphere grid using the original VCS profiles from iLemkd (119971 ) 
but with twice the value of the critical field. This grid is similar to that used by the Montreal 
group in the past ten years or so, and which has been applied to several studies of DA white 
dwarfs using the spectroscopic technique. 

Once again, we rely on the PG sample of DA stars analyzed above. A comparison 
similar to that shown in Figure [TT] is not very instructive here since the differences in T^s 
and logg are considerably smaller. Instead, we use the representation displayed in Figure \T2\ 
where the differences ATefj and A log g are shown as a function of effective temperature. This 
comparison reveals that our improved line profiles yield higher effective temperatures and 
surface gravities, with an important correlation with Tgfj. In particular, the logg values are 
about 0.05 dex larger above 20,000 K but can be as much as 0.1 dex larger near 15,000 K. 
Similarly, the differences in temperatures reach a maximum near 25,000 K but decrease 
at both lower and higher effective temperatures. The corresponding mass distributions, 
displayed in Figure [131 indicate that the mean mass is shifted by +0.034 Mq when our new 
models are used. The shape of the mass distributions is statistically equivalent, however, 
an d the dispersion remai n s the same. We must also point out that according to the results 
of iTremblay &: Bergeron! (120081 ). 15% of the DA white dwarfs in the temperature range 
considered here probably have thin hydrogen layers (Mn/Mtot < 10"*^). Consequently, the 
typical values for the mean mass are probably ~ 0.005 to 0.01 Mq lower than the values 
reported in Figure [T31 

As discussed in the Introduction, the spectroscopic technique provides very accurate 
measurements of the atmospheric parameters, allowing for a relative comparison of these 
parameters among individual stars. However, the absolute values of the atmospheric param- 
eters may suffer from an offset due to uncertainties in the physics included in the model 
calculations. Hence it is very important to compare the results with independent methods. 
The most reliable independent observ ational constra i nt for DA white dwarfs comes from 
trigonometric parallax me asurements. iHolberg et al.l (l2008al ) have shown using the recent 
photometric calibrations of iHolberg &: Bergeron! (120061 ) that there exists a very good correla- 
tion between spectroscopically based photometric distance estimates and those derived from 
trigonometric parallaxes. Here we compare absolute visual magnitudes instead of distances. 
We first combine trigonometric parall ax: measurements with ma gnitudes to derive Mv{tt) 
values. We then use the calibration of IHolberg fc Bergeron! (120061 ) to obtain My (spec) from 
spectroscopic meas urements of Tpff a nd log q. We selected 92 DA stars with known parallaxes 
from the sample of iBergeron et al.! (12007! ) (the uncertainties on the parallaxes must be less 
than 30%). The My values obtained from both model grids described in this section are 
compared in Figure [TH The agreement is very good within the parallax uncertainties for 
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both model grids. Hence, despite the fact that our new models yield higher values of Tefj 
and log (7 (i.e., smaller radii), these two effects almost cancel each o ther and the predicte d 
luminosities (or My) remain unchanged, and so are the conclusions of iHolberg et al.l (l2008al ). 
Another important constraint is provided by the bright white dwarf 40 Eri B for which a 
very precise trigonometric parallax and visual magnitude have been measured by Hipparcos. 
These measurements yield My = 11.01 ± 0.01, while we predict My = 11.02 ± 0.07 and 
10.97 ±0.07 based on our new and VCS models, respectively. Although both determinations 
agree within the uncertainties with the observed value, our new grid provides an exact match 
to the measured My value. 

We also point out that trigonometr ic parallax measurements a re available mostly fo r 
cool white dwarfs (iBergeron et al.l 119971 ) and nearby white dwarfs (IHolberg et al.l l2008bf). 



which are e xcluded from our anal ysis. In contrast with our approach here, iBergeron et al. 



(119971 ) and IHolberg et al.l (l2008bl ) make use of photometric measurements that cover the 
full spectral energy distributions, and both studies find mean masses near 0.65-0.66 Mq, 
consistent with our spectroscopic determinations for the PG sample. 



5. CONCLUSION 

We have combined for the first time in a consistent physical framework the unified 
theory of Stark broadening from Vidal, Cooper, & Smith and the HM88 non-ideal equation 
of state. Both of these well known theories represent the basis of our current knowledge of DA 
white dwarfs. Following the suggestion of the late Mike Seaton, we have taken into account 
the non-ideal effects due to proton and electron perturbations directly into the line profile 
calculations. We have shown that our improved line profiles are significantly different from 
other Stark profiles that have been used in previous studies. We have computed the first grid 
of model s pectr a without the need of the parameterization of the critical field introduced by 



Bergeron! (119931 ) and used in all previous models to mimic these non-ideal effects. We have 
demonstrated that our new profiles have important astrophysical implications. In particular, 
the mean mass of DA white dwarfs is ~ 0.03 Mq higher than previously measured. Yet, 
our updated atmospheric parameters determined from the spectroscopic technique remain in 
excellent agreement with the constraints imposed by trigonometric parallax measurements. 

Future work will confront our improved models with observations in the more complex 
regime of higher and lower effective temperatures. An exhaustive look at the Lyman line 
analysis of UV observations from FUSE is also much needed since the model spectra are 
highly sensitive to the non-ideal effects in that particular spectral region. For cooler white 
dwarfs, we must also investigate the abrupt cut-off that has been introduced to limit the 
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effects of the Lyman pseudo-continuum opacity. 
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Table 1. Opacity Sources Included 



Opacity 


Populations 


Cross Section 


H bound-bound, Stark broadening (H— H"*", H— e~) 


N«(H) 


This work 


H bound-bound, quasi-molecular (H-H, H— H+) 


N{H) 


AUard et al. (2004) 


H bound-bound, neutral broadening (H-H) 


N,(H) 


Ah & Griem fl965, 1966) 


H bound-free 


N,(H) 


Mihalas (1978) 


H free-free 


NeN(H+) 


Mihalas (1978) 


H2 free-free 


NeN(H+) 


we assume that of H I 


H3 free-free 


NeN(H+) 


we assume that of H I 


H~ bound-free 


N(H-) 


_John (1988) 


H~ free-free 


NeN(H) 


John (1988) 


HJ bound-free 


N(H+) 


Kurucz (1970) 


hJ free-free 


N(H)N(H+) 


Kurucz (1970) 


H^ free-free 


NeN(H2) 


Bell (1980) 


CIA H-H2 


N(H)N(H2) 


Gustafsson & Frommhold C2003) 


CIA H2-H2 


N(H2)N(H2) 


Borysow ct al. (2001) 


Rayleigh H 


N(H) 


Kissel (2000) 


Rayleigh H2 


N(H2) 


Dalgarno & Wilhams (1962) 


Thompson e scattering 


Ne 


Mihalas fl978) 



Note. — The bound-bound H2 and H3 opacities are always negligible in DA white dwarfs. 



- 29 - 



Fig. 1. — Theoretical line profiles of models at different effective temperatures and surface 
gravities using our new line profiles discussed in § 3. The lines correspond to H/5 to H8 
{bottom to top) from the Balmer series of the hydrogen atom. In each panel, the line profiles 
range from logg = 7.0 {dashed line) to 9.0 {thick line) by steps of 0.5 dex. The profiles have 
been convolved with a 6 A FWHM Gaussian profile, normalized to a continuum set to unity, 
and offset vertically from each other for clarity. 

Fig. 2. — Solutions in a T^g — logg diagram for a typical DA star using 1 line (H/?), 2, 3, 
4 and 5 lines (up to H8) in the fitting procedure (represented by thicker 1 a uncertainty 
ellipses from our fitting procedure). The top panel shows the results with the VCS line 



profiles, while the middle panel also includes the ad hoc parameter proposed by [Bergeron 



(I1993I ). The bottom panel is with our new line profiles discussed in this work (§3). 



Fig. 3. — Total occupation probability Wn for the bound states of hydrogen with principal 
quantum number n for T = 10,000 K, logNe = 17 {filled circles) and T = 20,000 K, 
logA^'e = 16 {open circles). 

Fig. 4. — Atomic hydrogen opacity at the photosphere of a 20,000 K DA white dwarf. The 
contributions of line and pseudo-continuum opacities are shown as solid and dotted lines, 
respectively. 

Fig. 5. — Stark broadening profiles as a function of wavelength measured from the line center 
at T = 10,000 K and logiVg = 17. Results are shown for H7 {left) and H8 {right). In each 
panel, we compare the results from this work (based on the VCS theory coupled with non- 
ideal effects from HM88; solid lines) with t he original VC S calculations {dashed lines). Also 



shown are the approximate calculations of ISeatonI (Il990l . dotted lines), which also include 
non-ideal effects. 

Fig. 6. — E r aissiy ity of a hydrogen plasma {filled dots) for the low-density experiment of 



Wiese et al.l (1l972l ). Our best fits {solid line) to the data are shown for different line broad- 
ening theories identified at the top of each panel along with the plasma state parameters 
obtained from the minimization procedure. 

Fig. 7. — Same as Figure 6 but with the data set for the high-density experiment of 



Wiese et al.l fll972l ). 



Fig. 8. — Corrections that must be applied to transform the atmospheric parameters obtained 
from the VCS profiles to our improved models. These have been obtained by simply fitting 
the Balmer line profiles of the VCS model grid with our new spectra. For clarity, we have 
omitted the grid points at logg = 7.75 and 8.25. 
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Fig. 9.— Left: For each star in the PG sample with 40,000 K > Teg > 13,000 K, we 
computed the average absolute deviation in T^s between solutions obtained from fits that 
include 2 to 5 lines (similar to Fig. [2]). These deviations were then averaged for all DA stars 
in our sample [filled dots). The results are shown for the three grids discussed in the text 
and identified on the x— axis. For comparison, the mean uncertainty of the fitting procedure 
is shown as dotted lines. Right: Same as left panel but for the dispersion in log g. 

Fig. 10. — Our best fit to the Balmer lines of WD 0205+250 with our improved line profiles 
[left panel) and with the VCS profiles [right panel). The atmospheric parameters are given 
in each panel. 

Fig. 11. — Mass versus T^s distribution for the DA stars in the PG sample in the range 
40,000 K > Teff > 13,000 K. Results are shown for both our improved lines profiles [top 
panel) and the VCS profiles [bottom panel). Lines of constant mass at 0.55 and 0.70 Mq are 
shown reference. 

Fig. 12. — Differences in Teg and logg obtained with our improved models and with the VCS 
profiles [Pcrit X 2) as a function of effective temperature for the DA stars in the PG sample 
in the range 40,000 K > Tcs > 13,000 K. The solid lines represent the 1:1 correlation while 
th e dotted lines c orresp ond to the uncertainties of the spectroscopic method as determined 



bv lLiebert et al.l (1200 



3SpO 



Fig. 13. — Mass distributions for the subsample of PG stars studied in Figure [T^ The mean 
masses are reported in the figure. 

Fig. 14. — Comparison of the absolute visual magnitudes obtained from trigonometric paral- 
lax measurements and from the spectroscopic technique using both model grids. The analysis 
is restricted to 40,000 K > Tes > 13,000 K. The error bars represent the ave rage uncertain 



ties. WD 1606+422 [open circle) is a suspected double degenerate binary (IBergeron et al. 



200l|). 
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Figure [H 
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Figure [2] 
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Figure [5] 
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